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ABSTRACT 

Core Accretion (CA), the de-facto accepted theory of planet formation, requires for¬ 
mation of massive solid cores as a prerequisite for assembly of gas giant planets. The 
observed metallicity correlations of exoplanets are puzzling in the context of CA. While 
gas giant planets are found preferentially around metal-rich host stars, planets smaller 
than Neptune orbit hosts with a wide range of metallicities. We propose an alternative 
interpretation of these observations in the framework of a recently developed planet 
formation hypothesis called Tidal Downsizing (TD). We perform population synthesis 
calculations based on TD, and find that the connection between the populations of 
the gas giant and the smaller solid-core dominated planets is non linear and not even 
monotonic. While gas giant planets formed in the simulations in the inner few AU 
region follow a strong positive correlation with the host star metallicity, the smaller 
planets do not. The simulated population of these smaller planets shows a shallow 
peak in their formation efficiency at around the Solar metallicity. This result is driven 
by the fact that at low metallicities the solid core’s growth is damped by the scarcity 
of metals, whereas at high metallicities the fragments within which the cores grow 
contract too quickly, cutting the core’s growth time window short. Finally, simulated 
giant gas planets do not show a strong host star metallicity preference at large separa¬ 
tions, which may explain why one of the best known directly imaged gas giant planet 
systems, HR 8799, is metal poor. 
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1 INTRODUCTION 

Given the be wildering diversity of exoplanet system arch i- 
tectures (e.g.. iBatalha et al]|2013l : IWinn fe Fabrvckvll2014h . 
it is clear that the outcome of planet formation is a highly 
stochastic process. Any statistical trends found in the ob¬ 
served planetary populations have special significance as 
they testament to planet-forming processes so robust that 
they rise above the stochasticity. A successful planet forma¬ 
tion theory must reproduce such trends. 

Metallicity correlations of the observed planets is one 
such correlation. Giant planets are detected much more fre- 
quent ly around metal-rich stars than arou nd metal-poor 
ones dGonzaledll999l : iFischer fe Valentill2005h . This correla¬ 
tion has been argued to provide a direct support to Core Ac - 
cretion theory for planet formation ('e.g. JPollack et ai1ll996l f 
since metal-rich en vironments assemble m a ssive cores much 
more readily (e.g., Ilda fe Linl 120041 . 120081 ; iMordasini et al.l 
12009k These cores are the crucial step towards forming 
a gas giant planet in CA framework, thus a positive gi¬ 
ant planet-metallicity correlation ensues. Gravi tation al disc 
Instability model for planet formation (e.g., I Bos ] 1 19981: 
iHelled et al.ll2013T ) was argued to produce a negative rather 


than positive giant planet-metallicity correlation since plan- 
ets’ radiative cooling is the fastest at low dust opacities 
(IHelled fe Bodenheimerll201lh . increasing planet’s chances 
of survival. 

However, there are both observational and theoretical 
reasons to study the issue further. Radial velocity observa¬ 
tions show that Neptune-mass planet occurrence does not 
seem to correlate with the host star’s metallicity for FGK 
stars (ISousa et al.ll2008l ). Similarly, transit method observa¬ 
tions with Keple r , sensitive to the planet’s radius rather than 
its mass, show (IBuchhave et al.ll2012l l that planets smaller 
than 4A® form around hosts with a wide range of metallic¬ 
ities, with the average close to the Solar metallicity. This is 
clearly surprising in the context of CA, since these smaller 
planets are the precursors of the gas giant planets. If CA’s 
explanation for the positive gas giant planet correlation with 
metallicity is correct then there should be more cores at 
higher metallicities, which is not what is observed. 

On the theoretical side, several important extensions 
to the classical variant of the Gravitational Instability (GI) 
model for planet formation have been recently proposed. 
In particular, GI model that includes planet migration and 
pebble accretion was showed to produce a strong posi- 
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tive and not negat ive correlation with the star’s metallicity 
(Nayakshin 2015a b) for coreless gas giant planets. The goal 
of this paper is to extend this recent work on giant planets 
with cores and also on core-dominated (smaller) planets and 
to compare the results with the observations. 

The key reason why GI theory may require a major 
overhaul is the realisation that GI gas fragments can also 
migrate from ~ 100 AU all the way into the inner disc 
(IBolev et alj|20ld l. While planet migr ation was a stan dard 
feature for the CA model since 1996 dLin et al.l Il996l l . GI 
planets were somehow thought to not migrate until recently. 
If this were true than GI could at best account for important 
but r are giant pla net sy stems like HR 8799 imaged directly 
fe.g.. iMarois et aLll2008r ). and could never explain the Solar 
System: massive self-gravitating proto-planetar y discs can 
hatch clumps only a t R Si many tens of AU (e.g. jRice et al.l 
120051: lRafikovll2005l'l . 

_ H oweve r, simulations, starting with IVorobvov fe B asu 

SB, showed t hat massive g as fra gments do migrate 
inward. Furthermore, IBolev et al.l (l 20 ldl suggested a new 
way of forming terrestrial-like planets inside ~ a few Jupiter 
mas ses gas clum p s by gr ain growth and s edimentation (see 
also IWilliams fe Crampinl[l97ll : [b osdl 19971 ), and then releas¬ 
ing them back into the disc by de stroying the gas c lumps 
via tidal forces from the host star. iNavakshinl (I2010al 'l used 
analytical estimates of the relevant processes and arrived at 
similar ideas, proposing the Tidal Downsizing (TD) hypoth¬ 
esis for formation of all types of planets observed so far. 

The key question for a planet forming in the framework 
of the TD hypothesis is how does the clump’s inward mi- 
;ration time scale, which may be as short as ~ 10 4 years 
Baruteau et al . 20 111 : iMichael et al.l l201ll ) , compare with 
the time scale for its internal evolution? If the fragment con¬ 
traction time is shorter than the inward migration time, then 
the fragment collapses, becoming a bona fide ’’hot start” 
protoplanet, which may mature into a gas giant planet. If 
the fragment contraction time is long, then the host’s tidal 
forces catch up with the fragment when the latter migrates 
too close to the star, and the fragment is disrupted. The 
fragment’s short existence may however not be all in vain as 
there could be a remnant. If grains within the fragment had 
a sufficient time to grow and sediment to the centre, getting 
locked into a self-gravitating massive core (see also KuijaeiJ 
1951; McCrea fe Williamsl 1965: Williams fe Crampin 197 it 


Boss Il997l : iHelled fe Schubertl 20081 1 . then the remnant of 


the disruption is the core - a practically ready rocky planet. 
No planetesimal accretion is required for the planet to sur¬ 
vive at t hat stage , althoug h ’ ’venee r accretion” of them or 
pebbles djohansen fe Lacerda|[201C)l ) is possible. 

Making detailed predictions in the TD hypothesis set¬ 
ting worthy of comparison to modern observations of ex¬ 
oplanets is however not trivial. Until recently, it was as¬ 
sumed that collapse of the self-gravitatin g GI gas fragments 


Bodenheimer et al j 198C 

: Cameron et al. 1982: Helled et al. 

20081: Navakshiii 201 ll: 

Navaksliin & Cha 20131'). The for- 


mation of giant gas planets in thi s pict ure is similar to 
that of low mass stars (e.g., lLarso S), except that ac¬ 
cretion of gas ont o the pre-c ollapse planets is terminated 
very early on (e.g.. lNavakshinl|2010bl ). This formation chan¬ 
nel f avors survival of mas sive (M p Si a few Mj) gas gi¬ 
ants dForgan fe Rice! 120131 ) , as more massive planets cool 


more rapidly (e.g., see lNavakshinl2015al 'l . Also, radiative col¬ 
lapse of ^-dominated fragments pr oduces a negative planet 
frequ ency - metallicity correlation (IHelled fc Bodenheimen 
dml) because radiative cooling is more rapid at l ow metal- 
licities. That is clearly at odd s with observations rtGonzaled 
Il999l : iFischer fc Valentill2005l ). 

We (INavakshinl l2015al ). however, considered accre- 
tion of 1 cm sized grains (usually called ’’pebbles”; 
l.lohansen fc Lacerdal 120101 : lOrmel fc Klahrl 1201(1) onto the 
pre-collapse fragments. While gas has pressure gradient 
forces that may prevent its accretion from the proto¬ 
plane tary disc onto the fragments (e.g., iNavakshin fe Cha 
12013) . pebbles weakly coupled to the gas by aerodynam¬ 
ical forces can still accrete onto the fragment in an anal- 
ogy to the pebble accretion pro cess onto massive solid cores 
dLambrechts fc Johansenl I2012T) in the context of the CA 
model. 

Surprisingly, pebble accretion on GI fragments was 
found to accelerate their contraction despite increasing their 
metallicity and hence opacity. Physically, addition of peb¬ 
bles to the fragment increases its weight without increasing 
its thermal energjQ. Since pre-collapse molecular gas frag¬ 
ments are polytropes with an effective polytropic index 7 
edging ever closer to the unstable value of 4/3 as the cen¬ 
tral temgeratureoUthe fragment increases towards 2000 K 
fsee lNavakshin|[2015al ). addition of a relatively small, ~ 10%, 
amount of mass in metals to the fragment usually tips it over 
into collapse. Higher metallici ty environments pr ovide larger 
pebble accretion rates, hence INavakshinl d2015bl) found, via 
a population synthesis like study of coreless gas fragments, 
that the frequency of gas giants survived the initial tidal 
disruption increases as metallicity of the host in creases. 

An important issue not considered by INavakshinl 
d2015al lbh is formation of a core within the fragment. If peb¬ 
bles entering the fragment accrete onto the core, the core’s 
accretion luminosity may be significant and this may lead 
to expansion of the fragment instead of its contraction. A 
similar effect - accretion of planetesimals onto the grow¬ 
ing core - is well known in the context of the CA model 
and has to be take n into account when determining the fat e 
of the planet (e.g. , IPollack et al .1 117)961 : lAlibert et alj|2005h . 
INavakshinl d2015d . paper I hereafter) has recentl y presented 
numeri cal algorithms that extend the work of iNavakshiid 
d2015bl ) by adding the core formation and growth inside the 
pre-collapse fragments. 

Here we use that framework to study the planet fre¬ 
quency of survival versus metallicity dependence for all types 
of planets rather than just the coreless giants. We also go 
beyond of our previous work by studying not just the out¬ 
come of the migration vs disruption competition, but also 
how the population of planets formed in TD model would 
look like in the planet mass - separation plane. 

The structure of our paper is as following. In (J2] we 
present initial conditions, assumptions and numerical algo¬ 
rithms that are employ to calculate an outcome of a gas-dust 


1 This is in difference to accretion of planetesimals in the CA 
framework, in which planetesimals impact the growing planet at 
a high (up to a few km s —4 ) velocity and hence may actually heat 
it. In contrast, pebbles sediment onto the fragment at a moderate 
~ a few m/sec velocity, or else they fragment. Their kinetic energy 
input into the planet is hence negligible. 
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Planet-metallicity correlations 3 


fragment being born in an outer region of a massive proto¬ 
planetary disc. In fj3l two example planet formation tracks 
are shown, one for a gas giant and another for a hot Super- 
Earth planet. In fj4] we present assumptions and methods 
for calculating a grid of models in a population synthesis 
like study. (0 presents the population synthesis study re¬ 
sult in the planet mass versus planet-host separation plane. 
In |j 6 ]we marginalise over the results to derive the planet- 
metallicity correlations for our simulated planets. Discussion 
of main ideas and results of the paper is given in m whereas 
conclusions are presented in fj 8 ] 


Initial disc properties, M d =0.15 M sun 



2 NUMERICAL METHODS 
2.1 Disc treatment 


We follow the procedures described in paper I with several 
changes detailed now. Our goal is to simulate formation of 
planets from their birth to the end of the disc migration 
phase. The initial disc surface density profile is given by 



where Ro is the disc length-scale, set to Ro = 100 AU, 
Rin = 0.08 AU is the inner boundary radius. The constant 
Am is calculated so that the disc contains a given initial 
mass Md = 2% J^ out RdR'EoiR)- The disc is ’’live”, that 
is evolved by solving the time-dependent ID viscous disc 
evolution equations that include the planet-disc interactions 
(section 3 in paper I), and now also include the disc photo¬ 
evaporation as a sum of the UV and the X-ray driven terms, 
using the fits to the p h oto-ev apor ation rates S e v(R) from 
lAlexander fe Armitagd d2007t) and lOwen et al.l d^OlT ). re¬ 
spectively. The ionising photon luminosity of the star is set 
to $ion = 10 42 photons s -1 , while the X-ray flux of the star 
is Lx = 2 x 10 30 erg s -1 . 

Gas fragments in the disc are born by self-gravitational 
instabilities, hence the disc must be initially gravitation¬ 
ally unstable in its outer region. Figure [T] shows the ini¬ 
tial disc properties for disc mass Md = 0.15 Mq and stel¬ 
lar mass M» = IMq. The top panel shows the disc sur¬ 
face density, S(R), as well as the disc photo-evaporation 
rate profile (which does not change during the simulation). 
The quantity plotted with the red dotted curve is X e v7r R 2 , 
in units of 10 -11 Mq year -1 . The bottom panel shows the 
Toomre Q-parameter (black solid curve) and the dimension¬ 
less cooling time, t* ool = t coo iQ(R). For the disc to be self- 
gravitating and actually fragmenting, Q ^ 1.5 and t* ool ^ 3 
dGammiel 120011 1. Both of these conditions are satisfied at 
R ~ 70 — 80 AU. 

The disc v isc osity is described with the 
IShakura fe Sunvaevl d 19731 ') viscosity parameter, but 
we now add a self-gravity contribution to it: 

(X — CXq -f- Cksg , (2) 


where Qo is radius and time-independent (free) parameter 
of the model, whereas a a g describes gravito-turbulence 


CCag 


= 0.2 


Qo 

Ql + Q 2 ' 


(3) 



Figure 1. Initial disc properties for Md = 0.15 Mq. Top: Ef R) 
and the wind mass loss rate profile, defined as Y,(R)nR 2 , as la¬ 
belled. Bottom: Toomre parameter, Q, and dimensionless cooling 
time of the disc, t coo lfl(R), versus radius. The disc is gravitation¬ 
ally unstable and may fragment at R ~ 70 — 80 AU. 


where Q o = 2 and Q is the local Toom re’s parameter. This 
form is inspired bv lLin fe Pringld dl987l ~l suggestion that a oc 
Q~ 2 in the self-gravitating regime, on the one hand, and 
mor e recent simulati ons showing that a sg saturates at about 
0.1 dRice et al.ll2005l ~). on the other. 


2.2 Planet’s birth and migration 

The following considerations are important for both the ini¬ 
tial and the final phases of our simulations. Simulations of 
self-gravitating protoplanetary discs show that fragments 
rarely form in isolation, hence we believe that every star 
is likely to have hatched a number of fragments early on. 
Presumably, most of these fragments migr ate rap idly and 


Presumably, most ot these fragments migr ate rap idly and 
perish by being disrupted (see, in particular. I Vorobvovl 20131 ') 
or by being driven all the way into the star. This picture is 
consistent with the suggestion that disruption and swallow¬ 
ing of planets by thei r host stars power FU Ori outbursts 
of youn g protostars (IVorobvov fe Basul 120061 : iBolev et all 
l2010l : iNavakshin fe Lodatol 20121 '). statistics of which sug- 
gests that every star may have ~ 10 — 20 of such episodes 
dHartmann fe Kenvon|[l996ll . 

Given that the o bserved giant gas planets ’ frequency of 
occurrence is 10% dWinn fe Fabrvckvll2014lf per star, the 
survival of a giant gas planet is a very rare event indeed. It 
is reasonable to assume that giant gas planets that survive 
to be observed are those that were either hatched by the 
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4 S. Nayakshin 


disc late, or those that hang around at the outer disc for a 
while before migrating in, because such fragments are more 
likely to avoid the tidal disruption and also being driven all 
the way into the star. 

This view is not unreasonable. Simulations show that 
presence of one clump in a disc may induce for mation 
of more clumps in the same disc (e.g., iMer ul l2013ll . Mul¬ 
tiple clump formation leads to strong interactions be¬ 
tween the clumps, so that some of them migrate inward 
faster while others ar e scattered on larger orbits (e.g., 
ICha fe Navakshinll201ll ): some of the clumps may be com¬ 
pletely ejected from the system by the clump-clump inter¬ 
actions ( Basu fe Vorobvovll2012l !. Clumps on scattered but 
still bound orbits may take a while to loose their orbital 
inclinations with respect to the midplane of the disc, and 
eccentricities, and only then start migrating in. 

Furthermore, presence of clumps in the disc may in 
fact change the fragmentation properties of the disc since 
the c lump induces additional perturbations (see again [Merul 
boil , where discs were found to fragment at surprisingly 
small radii i n the presence of a pre- existent gas clump at 
larger radii'l. iNavakshin &: Chal (120131 ') focused on migration 
of a single fragment in a marginally self-gravitationally sta¬ 
ble disc, thus placing the fragment in a massive but stable 
disc that would not fragment on its own. They found that 
in the presence of the clump the disc became more unstable 
and formed additional fragments. This shows that it is pos¬ 
sible to hatch new gas clumps in discs with Q > 1.5 if there 
are already fragments born earlier. 

Clearly, when survival of the fragments is rare, we 
should not under-estimate the importance of initial condi¬ 
tions giving the clumps the best chance of surviving. There¬ 
fore, our population synthesis model aims to account for this 
by (a) slowing down the type I migration below its nominal 
rate by some factor larger than unity (see below); (b) allow¬ 
ing the disc mass to be below the Md = 0.15 Mq marginally 
gravitationally unstable configuration when the fragment is 
born, and (c) removing the disc instantaneously after a time 
t r em randomly sampled between 0.5 and 5 Million years. 
These choices are certainly not unique but we feel they are 
reasonable. An extremely rapid external phot o-evap oration 
of the disc by a close-by massive star dClar kel bOOTt) or the 
presence of a secondary at beyond 100 AU that is later re¬ 
moved by star-star interactions may be important for such 
a rare population of planets as the gas giants. 

The disc exchanges angular momentum with the planet 
via type I and/or type II migration torques. The only dif¬ 
ference from paper I is that the type I migration time scale, 
ti, is now given by 


tl — /mig 


Mi 


MpMd H 2 


a n _1 


( 4 ) 


The factor (1 + M p /Md), used in paper I on the right hand 
side of this equation, has been dropped, as it is argued 
that the saturation of the type I migration rate in the limit 
M p /Md 1 is automatically taken into account by passing 
the torque of the planet to the ID viscous disc evolution 
code. We find that when this torque is large then a gap may 
start opening, thus reducing the torque. 

In equation @ /migr > 1 is a free parameter that ac¬ 
counts for a slower migration rate of the planet in a non self- 
gravitating disc compared to the isothermal type I migration 


rate. Since our discs are sampling the post self-gravitating 
phase of the disc evolution by design (see below), we may ex- 
pect type I migration to be slower than that for Q ~ 1 discs 
dBaruteau et al.ll201ll l. for which /migr ~ 1. In the models 
below /migr is randomly sampled between 1 and 10. 


2.3 Planet’s internal evolution 

Planet’s internal evolution is calculated as in paper I, util¬ 
ising the "follow the adiabats” approximation to convec¬ 
tive/radiative cooling of the planet. The planet is initialised 
as a polytropic sphere of gas of homogeneous composition, 
with metallicity, zo, equal to that of the parent disc, with 
a given central temperature To (see ©■ The planet cools 
b y emission of r adiation assuming interstellar dust opacity 
dZhu et al.ll2009ll multiplied by the factor / op < 1 to account 
for grain growth and also by the metallicity of the planet in 
Solar units, that is, by z/Z q. Since the planet is surrounded 
by the disc, its irradiation by the disc (or by the star if the 
planet is directly exposed to it in the final stages of calcu¬ 
lation) reduces the ra te at which the planet cools, creating 
a the rmal bath effect dCameron et al.ll 19821 : IVazan fe Helledl 

l2012|j . 

The planet accretes peb bles from the surrounding dis c 
at the Hill’s accretion rate dLambrechts &: Johansenll2012r ) 

M z =2f p E g (a)n a R 2 H (5) 

where Rh is the planet’s Hills radius, = (GM t / a 3 ) 1 ^ 2 , 
and £ g = / g E(a) is the grain surface density at radius R = 
a. The pebble mass fraction, 0 < / p < 1, is a free parameter 
independent of the planet’s location, fixed for a given run 
(see [|4] below). 

Note that pebble accretion rate depends directly on the 
surface density of gas around the planet. If the planet man¬ 
ages to open a deep gas in the disc, £(a) drops to very low 
values and hence M z dives to negligible levels. Pebble accre¬ 
tion is also turned off if and when the fragment undergoes 
the second (Hydrogen molecule dissociation) collapse. It is 
argued that after the collapse the planet is so compact that 
it should be able to accrete gas from the disc as well, so 
that both gas and pebbles would be accreted. We currently 
do not include gas accretion onto post collapse planets. 

Pebbles entering the fragment have a fixed initial size 
of a pe b = 0.03 cm for this paper, and are deposited in the 
outer several mass grids of the planet, where they are mixed 
homogeneously with the grains already in those regions. 
The grains are allowed to grow by sticking collisions and 
sediment into the centre of the fragment into the core, as 
described in paper I. Turbulence within the fragment, pa- 
rameterised by a coefficient ctd, and convection, however, 
tend to counteract grain sedimentation. We usually find 
that grain size must increase to about a g ~ 1 cm before 
they sediment efficiently. A further significant obstacle to 
a rapid grain sedimentation into the core is the fact that 
grains fragment in collisions rather than stick if collision 
velocity is too high, as shown by laborat o ry experiments 
on dust growth (e.g..lBlum fc Munchlll993l : iBlum fe Wurml 
120081 : iBeitz et al.ll201lf) . As explained in paper I, section 4.3, 
grain growth turns into grain fragmentation if grain sedi¬ 
mentation velocity with respect to the surrounding gas ex¬ 
ceed Umax, set in this paper to 3 m s -1 . Grains are also 
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vaporised if the surrouding gas temperature is too high for 
the given grain species. 

Three grain species are included in the model: wa¬ 
ter ice, CHON, and silicates ("rocks”). The relative abun¬ 
dances of the species are 0.5, 0.25 and 0.25, respectively. The 
grain growth, fragmentation, vaporisation, and sedimenta¬ 
tion equations are solved for each species separately since 
these three species have very different material and thermal 
properties. 

Grains reaching the centre are allowed to accrete onto 
the solid core there, whose initial mass is set to a very 
"small” value ( 1 O - 4 M 0 ). Growing core is expected to 
radiate some of its gravitational potential energy away, 
but a self-concistent modelling of energy transfer within 
the core is fraught with many physical uncertainties (e.g., 
iStamenkovic et all 120121 '). As in paper I (section 4.4), we 
parameterise the energy release by the core via the Kelvin- 
Helmholtz contraction time, tkh, of the solid core, which is 
fixed for all the runs here at f m in = 3 x 10 5 years. The lu¬ 
minosity released by the core is an internal energy source 
for the contracting fragment, analogous to that of a stel¬ 
lar core, and may even exceed the luminosity of the whole 
planet, causing its expansion. 


3 EXAMPLE RUNS 

Before the grid of models is discussed, we present two ex¬ 
ample planet formation tracks that illustrate the sequence 
of events in the TD model for planet formation. 

3.1 A gas giant planet 

Figure [2] shows an example selected from the grid of runs 
in which the end result is a gas giant planet of mass M p = 
1.17 Mj at separation of a = 2.7 AU. The initial mass of the 
planet is M p = 1.0Mj, so all the extra mass in the end of 
the run comes from accretion of pebbles from the disc. The 
initial mass of the disc is Md = 0.074 Mq, planet location 
ao = 83 AU, and the initial planet’s central temperature 
T c = 254 K. The disc metallicity is Solar, that is, 2 = 0.015, 
the pebble mass fraction is / p = 0.086, disc viscosity pa¬ 
rameter ao = 2.9 x 1CU 3 . The disc removal time is t rem = 3 
Million years, but in this simulation the disc dissipates ear¬ 
lier due to photo-evaporation, as we shall see below. Further, 
the opacity reduction factor / op = 0.5, the planet’s type I 
migration speed is moderated by factor /migr = 8.05, which 
is on a high side for the grid of models, so this run presents 
an example of a rather slowly migrating clump. 

Fig. EH presents the disc surface density profile at five 
different times. The initial condition at t = 0 is shown with 
the black solid curve, all the other curves are marked by their 
respective times in the legend. The crosses on the bottom 
of the panel show the position of the planet at times colour- 
coded in the same way as the surface density curves. 

Next panel, fig. (2jj, reports the radial location of the 
planet (black solid curve), the planet’s Hill (red dotted) and 
its actual (R p , blue dashed) radii, all ploted versus time. 
Boxed text and the corresponding arrows mark the times of 
three important transitions in this run. In particular, the 
fragment is initially molecular Hydrogen dominated, col¬ 
lapses at time t = 0.55 Million years, at which point Hydro- 


Disc evolution 




Figure 2. Evolution of the disc (panel a) and the embedded 
fragment (panles b and c) that survives to be a gas giant planet 
(see 13. 11 lor detail). Panel (a) shows disc surface density profiles 
at times t = 0 (solid curve) plus several later times as labelled in 
the legend. The position of the planet at corresponding times is 
marked by a cross of same collor at the bottom of the panel. Panel 
(b) shows the planet’s separation, radius and the Hills radius, 
whereas panel (c) shows the mass of the core versus time. 

gen becomes atomic and partially ionized. During collapse 
the fragment contracts to the radius of just a few times that 
of Jupiter, and then continues to contract further due to a 
rapid radiative cooling. This fragment is never challenged 
by tidal forces of the star as the planet’s Hill radius (red 
dotted curve) is always much larger than R p . 

The next notable transition in the planet-disc system 
occurs at t ft 1 Million years, when the planet opens a deep 
gap in the disc. The planet is located at a ft: 27 AU at 
that time. From that time on, the planet migrates in type 
II regime. Note that at time t = 1.32 Million years (the 
dashed blue curve in fig. [2H), the gap becomes deeper than 
it was before, and the inner disc has a significant gas deficit 
at a ~ 2 AU. This is not only due to the gap openning; the 
photo-evaporation of the disc is especially strong at a ~ a 
few AU, and is negligible within 1 AU, as can be seen from 
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fig. 1. The planet-facing side of the disc inside the planet’s 
orbit thus loose mass rapidly due to photo-evaporation. 

At this point in time, the inner disc is drained roughly 
equally efficiently by both accretion onto the star and photo¬ 
evaporation. By time t = 1.5 Million years, there is a com¬ 
plete gap between the planet and the star. The outer disc is 
then evaporated too, and the planet stops migrating com¬ 
pletely. 

Finally, fig.^Jj shows the core mass assembled within the 
fragment, and its decomposition on the silicates (“rocks”) 
plus water ice and CHON. Since the fragment heats up 
quickly when pebble accretion commences, water ice grains 
vaporise before they could sediment into the core, hence 
they do not contribute to the core at all. CHON grains 
(green dashed curve) sediment while the centre of the frag¬ 
ment is cooler than T ~ 400 K, and then only the rocks 
are refractory enough to continue to accrete onto the core. 
This particular planet did not manage to assemble a mas¬ 
sive solid core before it collapses (when all grains not yet 
locked into the core are vaporised); the final core mass is 
only M c = 0.55 M®. 

The final metallicity of the planet is z = 0.158, about 
ten times that of the initial (Solar) disc metallicity. This 
example thus shows how a metal-rich gas giant planet with 
a low mass core can be assembled in the context of the TD 
hypothesis. 

3.2 A hot Super-Earth planet 

Figure [3] demonstrates how a hot Super-Earth planet can be 
assembled in the framework of TD. The initial mass of the 
planet is again M p = 1.0 Mj. The initial mass of the disc is 
Md = 0.066 M®, planet’s location is ao = 108 AU, and the 
initial planet’s central temperature T c = 282 K. The disc 
metallicity is 0 = 1.7 times Solar, the pebble mass fraction 
is / p = 0.08, disc viscosity parameter ao = 2.3 x 10~ 3 . 
The disc removal time is t rem ~ 5 Million years, but again 
disc dissipates earlier due to photo-evaporation. The planet’s 
type I migration speed is barely reduced from the isothermal 
result, e.g., factor / m i gr = 1-3, so the fragment migrates 
in much sooner than it does in T3.ll The maximum grain 
sedimentation velocity for this simulation is set to u max = 
10 ms” 1 , which, combined with a higher 2 value, encourages 
a faster solid core growth. 

It appears that the much more rapid migration is to 
blame for most of the differences in the fate of the fragment 
here as compared to the much more slowly migrating one in 
H3.ll Fig. 3 shows that the fragment finds itself in the inner 
disc much quicker than the one in Fig. 2 did. The gap in 
the disc is openned by the fragment at around a = 7.2 AU, 
at time t = 0.18 Million years. Significantly, the fragment is 
still in the pre-collapsestate at that moment, and it also has 
a massive solid core, M c = 4.8 M® at the time of the gap 
openning. The more massive core is due to the higher value 
Umax in this simulation. 

Contraction of the fragment in this simulation occurs 
almost entirely due to pebble accretion because radiative 
cooling is very inefficient. This is because the fragment’s 
dust opacity is high, especially so due to pebble accretion: by 
t = 0.18 Million years the fragment metallicity is a = 0.144. 
The fragment is also basked in the radiation field of the disc 
(which at the outer disc radii is dominated by irradiation 


Disc evolution 




Figure 3. Same as Fig. [2] but for a simulation that forms a 
hot Super-Earth planet as the end result (E2). The blue box 
in the top panel (a) shows where and when the gas fragment is 
tidally disrupted, which leaves behind a massive M c ~ 6.4 M® 
core that continues to migrate, and finally arrives in the “hot 
region”, a = 0.23, by the time the disc is dispersed. 

from the star). Now, when pebble accretion onto the frag¬ 
ment stops, its contraction stops as well. However, since the 
fragment is metal-rich, the core of the fragment continues 
to grow in mass rapidly and its luminosity becomes suffi¬ 
ciently high to start puffing the fragment up. This increase 
in the fragment’s radis after the disc gap openning is easy to 
spot in Fig. 3b. Since the fragment continues to migrate in, 
now in type II regime, the planet’s Hills radius continues to 
shrink, whilst the planet’s radius, R p , increases with time. 
A tidal disruption is thus unavoidable, and indeed occurs 
when the fragment reaches a = 3.16 AU. 

The core’s mass is then equal to M c = 6.4 M®. As the 
fragment is disrupted, its mass plumets to M c (in this paper 
we neglect a possible tightly bound dense gas “atmosphere” 
around the core, which is not likely to be massive as our 
cores are quite bright). With this reduced mass, the planet 
is no longer able to maintain the deep gap in the disc, which 
is swiftly closed (see the blue dashed curve in Fig. |3ji). The 
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mass of the core is however sufficiently high for it to continue 
its migration in type I now, on a longer time scale yet fast 
enough to arrive at a = 0.23 AU by the time the disc is 
dispersed at t = 1.5 Million years. 

The composition of the core is almost exclusively rock, 
with CHON and especially water ice grains unable to sedi¬ 
ment because the core heats up to temperatures above ~ 400 
K quickly (cf. paper I for more detail). The core’s compo¬ 
sition is significantly different from that predicted by the 
Core Accretion paradigm, in which ices ar e very impor¬ 
tant d Pollack et al.l 1 19961 ; lAlibert et al.l l2005l l for assembly 
of cores as massive as this one. This may be a testable dif¬ 
ference between the models if composition of Super-Earths 
is constrained observationally, although the potential pres¬ 
ence of a Hydrogen/He atmosphere on top of the core may 
complicate differentiation between the models. 


4 THE POPULATION SYNTHESIS GRID OF 

MODELS 

A population synthesis involves investigation of the model’s 
parameter sp ace, which is usually large, by randomly sam¬ 
pling it fe.g.. IMordasini et al.ll2009r ) in order to pull out sta¬ 
tistical trends of the model that may not be distinguisheable 
from just a handful of runs. In following this approach, we 
Ex some parameters of the model in order to reduce the pa¬ 
rameter space to a manageable but physically meaningful 
minimum. 

Table 1 lists these fixed parameters, their meanings and 
values. The latter are chosen to be “reasonable” based on 
experiment (e.g., Umax) or on what is commongly used in 
literature (e.g., Lx)- While these values are not unique, we 
found by varying them that none of the parameters in Table 
1 change the metallicity trends that are the focus of this pa¬ 
per significantly, although quantitatively there are of course 
changes in the results. For example, the normalisation of 
the giant planet number per star surviving in the end of the 
simulations does depend on how the disc is removed in the 
end of the simulations, and hence on Lx and <3>i on , but the 
metallicity trends of the survived planets are the same. 

One parameter in Table 1 that is difficult to constrain 
from first principles and which does influences the results is 
fop, the grain opacity reduction due to grain growth. The 
latter is fixed here on a rather large value, / op = 0.5, which 
essentially precludes a significant opacity reduction. This 
is done following the arguments presented in paper I: the 
observed metallicity trends are incompatible with low dust 
opacity models and presumably imply that dust fragmen¬ 
tation (and not only dust growth) is occuring within the 
fragments to keep the opacity at relatively large values. 

Table 2 lists the set of variable parameters, and their 
minimum and maximum values. The parameters for a run 
are then randomly picked in the logarithmic space between 
these respective minimum and maximum values. For exam¬ 
ple, the disc viscosity parameter, ao, is very poorly known 
for protoplanetary discs (and any other astrophysical discs). 
We introduce a random variable £ a uniformly distributed 
between 0 and 1 , draw a random value and define ao for 
a given simulation by 

In a 0 = In a min + (1 - £«) In a max . ( 6 ) 


Table 1. Fixed parameters of the simulations and their values: 
maximum velocity for sticking rather than fragmentation, u ma x, 
in m s — , core’s Kelvin-Helmholtz contraction time, t^h, in years; 
opacity reduction factor due to dust growth, / op ; pebble radius, 
a P ebi in cm; ionising photon flux from the star, 4>; on , photons 
s — , X-ray luminosity of the star, Lx, in erg s —1 . 


Parameter 

%ax 

*kh 

/op 

^peb 

<f>. 

^ion 

Lx 

Value 

3 

3 X 10 5 

0.5 

0.03 

10 42 

2 x 10 32 


Table 2. Randomly varied parameters of the grid of models. The 
first row gives parameter names, the next two their minimum 
and maximum values. The columns are: ao, the disc viscosity pa¬ 
rameter; / p , the pebble mass fraction; factor; ao [AU], the initial 
position of the fragment; M<j, the initial mass of the disc, in units 
of Mg; /migr, the reduction multiplier for type I planet migra¬ 
tion; hem, Million years, the time for disc removal; ad, turbulence 
parameter within the fragment. 

Parameter a f p a 0 M A / m i gr hem ad 

Min 0.001 0.05 60 0.05 1 0.5 10“ 4 

Max 0.01 0.2 120 0.2 10 5 0.003 


Similarly, random uniformly distributed variables are gener¬ 
ated for all the other entries in Table 2. 

Having defined the parameters for a simulation, we 
model one gas fragment per disc at a time. The fragment 
starts with a given initial total (H/He + grains) mass, and 
has the same metallicity as that of the disc. The initial 
temperature is randomly picked (as described in equation 
©above), between To(M p ), given by 


Tq(M p ) = 100 K 


M p l 1 / 2 
0.5 Mj. 


(7) 


and twice To(M p ). Here the mass dependence of To reflects 
the fact that more massive fragments c ool more rapidly (e.g., 
I Helled et alj l200Sl : iNavakshlnl l2010bl ~l , and are also hotter 
at formation due to the adiabatic compression. The choice 
made in equation [7] does not pre-determines or drives the 
main results of this paper. 

Finally, we wish to investigate how our results depend 
on the metallicity of the host star/disc and the initial mass 
of the fragment. Therefore, the random models are repeated 
for a grid of the metallicity and the initial fragment’s mass 
values. The metallicity grid is given by Zi/z © = 3 b -3 )/ 2 , 
with i = 1,2, ... 5. The initial planet mass grid is given by 
M p / Mj = 0.5, 0.7, 1 and 2. For each of the combinations 
of z and M p from this grid, e.g., 2 = 3 Zq and M p = lMj, 
the total of 243 models are run. This gives the total of 4860 
runs. A typical run takes about 2 hours of physical time on 
a single CPU, although some runs (usually with higher disc 
mass and viscosity, as that leads to smaller time steps) had 
to be executed for up to 48 hours. 


5 RESULTS 

Figure [4] presents the results of the 4860 runs in the planet 
mass versus planet-star separation plane. In order to facil¬ 
itate a visual comparison to the exoplanetary data, giant 
planet masses are shown as M p sin i, where i is inclination 
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of the system to the observer, generated randomly assuming 
a uniform distribution for cosi. This is done because many 
of the known giant planets were detected by the radial veloc¬ 
ity method and so also have the unknown sin* factor. The 
gas giants are defined here as fragments that did not have 
a tidal disruption and hence are all more massive than our 
minimum starting fragment mass, e.g., 0.5 Mj. All of these 
planets are above the horizontal dashed line. 

Cores, the remnants of tidal disruption of the initial 
gas fragments, are the planets below the dashed line. In Fig. 
[2 all of the survived giant planets are shown, but, for the 
sake of clarity, only a third of the cores, randomly selected 
from the full sample, is plotted. The mass of the cores, the 
scale for which is given in units of M® on the right vertical 
axis, is not multiplied by the sini factor since many of the 
observed planets of this mass are detected by the transit 
method which is not affected by the inclination uncertainty. 

We emphasise that in this paper no consideration for 
the presence of a bound atmosphere around the cores af¬ 
ter the disruption is made, so their mass may be somewhat 
under-estimated. Further, gas accretion from the surround¬ 
ing protoplanetary disc onto the cores and onto the post¬ 
collapse giant gas planets may well occur but is also not 
taken into account, thus again possibly under-estimating the 
final mass of the planets. Finally, an additional uncertainty 
which may bring down the mass of the giant gas planets is 
that here we assumed that H 2 dissociation collapse of a gas 
giant leads to the instantaneous collapse of the whole frag¬ 
ment. However, 3D hydrodynamical simulations of embed¬ 
ded clumps in protoplanetary discs show that only a fraction 
of the total mass of the fragment collapses into the post¬ 
collapse planet imm ediately, with the res t (up to ~ 50% of 
the fragment’s mass iGalvagni et al.l 1201211 collapsing into a 
circum- planet ary disc. It is not clear whether all of this disc 
or only a fraction is eventually accreted onto the planet. 
There is thus an uncertainty of a factor of ~ 2 in the mass 
of the giant planets shown in Fig. [2 

Despite these uncertainties, there are clear and strong 
metallicity trends in the planet populations in Fig. [2 which 
we shall investigate in detail shortly. To aid this, the symbols 
in Fig.[2are colour-coded to show the metallicity of the host 
star/disc. For example, the red colour shows the most metal 
rich systems, 2 = 32®, whereas the green show the most 
metal-poor systems in the sample, z = (1/3)2®. It is easy to 
notice that most of the giant planets in the inner few AU of 
the diagram are metal rich; this is not so for cores, however. 
This observation essentially summarises the main findings 
of this paper. 

Finally, note the vertical column of giant planets at 
a ~ 0.09 AU. These are the fragments that managed to 
collapse before being tidally disrupted but that continued 
to migrate inward rapidly and reached the inner boundary 
of our computational domain (we stop the calculation if the 
planet reaches a short distance away from R[ n = 0.08 AU). 
This column contains 2391 planet, e.g., almost a half of all 
fragments born in the disc, managed to collapse before be¬ 
ing tidally disrupted but reached the innermost disc radius. 
Most of these would probably be completely devoured by 
the star if our computational domain were extended all the 
way to the inner boundary. 

An additional point to make on this is that our calcu¬ 
lation is optimised for survival of the giant planets because 


we essentially simulate the end phase of the disc evolution. 
There is no reason why fragments could not be born be¬ 
fore the phase we are simulating. These would migrate even 
faster, and most likely would all end up either being tidally 
disrupted or being driven all the way to the inner bound¬ 
ary. Hence the number of giant planets migrating all the 
way into the star is still an under-estimate. This is consis¬ 
tent with the view presented earlier - that destroyed gas 
giant planets cause FU Ori outbursts and that statistically 
speaking, there may be many such events per protostar. 


6 METALLICITY CORRELATIONS 


6.1 Small separations 

Figure [5] shows the number of planets formed in our runs at 
separations 0.1 AU < a < 5 AU (this excludes the planets 
that reached the innermost disc radius, e.g., the column of 
planets on the left-hand side of Fig. 0 versus host metallicity 
in Solar units, 2 / 2 ®, for three groups of planets: the gas 
giants, the Super-Earths and the “Earths”. The division on 
the three groups is made based purely on the planet’s mass. 
Gas giants are fragments that collapsed rather than were 
disrupted, so they are all more massive than 0.5 Mj. Super- 
Earths are solid cores more massive than 5M®, and Earths 
are planets with mass 1 M® < M v < 3 M®. 

There is a clear and strong correlation with 2 for gas 
gi ant planets , althou gh somewhat less strong than found 
in iNavakshin (|2015b|). As exp lained in the latter publica- 
tion and in iNavakshinl (l2015al l. physically this correlation 
arises from gas giants collapsing faster at higher abundance 
of pebbles in the outer disc, since the fragments then accrete 
pebbles at higher rates. A larger fraction of gas clumps thus 
collapses before it is tidally disrupted. 

The positive correlation with 2 in Fig. 0 saturates at 
the highest metallicity bin. It is worth emphasising that we 
assumed here that the fragment birth locations do not de¬ 
pend on metallicity of the disc for simplicity. Higher 2 mean 
higher dust opacity, and hence longer cooling times in the 
disc. One may expect the clumps to be born farther away on 
average at higher 2 due to higher disc opacity (such effect 
is actually found in simulations bv lMeru fe Batal20ld f. The 
farther away the clumps are born, the longer it takes for 
them to migrate in, and hence the more likely they are to 
survive. This secondary effect, not included in the present 
simulations, may strengthen the giant planet frequency ver¬ 
sus metallicity trend of the model. 

On the other hand, Fig. 0 shows, surprisingly at the 
first glance, no clear correlation between the host metallic¬ 
ity and the populations of Earths or Super-Earths. Naively 
speaking, one may expect that a higher abundance of met¬ 
als within the fragments at higher metallicities would cause 
more massive cores to be built. When some of the fragments 
are disrupted, these cores should become visible to the ob¬ 
server, so this would predict more massive cores around stars 
of higher metallicity. 

This simple argument does not actually work because 
the fragment’s properties also vary with metallicity and that 
should be also taken into account. In particular, a fragment 
loaded with pebbles at a higher rate indeed has a higher 
metallicity, so that the core’s accretion rate is higher, but 
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Figure 4. Planet mass versus planet-host separation for the pebble accretion grid of models for different metallicities, as shown in the 
legend in Solar metallicity units. For the gas giant planets (historically found by RV methods) M p sin i is used as proxy for M p , where 
sin i is randomly generated orbital inclination with respect to the line of sight. For core-dominated planets, located below the dash-dotted 
line, M p in units of Earth mass (on the right vertical axis) is plotted. 


the time window for grains to sediment down into the core 
is shorter because the fragment heats up too quickly and the 
grains are eventually vaporised which clearly stops accretion 
of the core. 

This physics is internal to the fragment and has noth¬ 
ing to do with the fragment-disc interaction. Therefore, to 
investigate it in greater detail, we perform a series of runs in 
which we turn off fragment’s migration and keep the frag- 
ment’s pebble acc retion rate constant, parameterised, as in 
iNavakshinl (l2015al ). in terms of the ” metal loading time”, t z , 
defined by 

dM z zqM p , 

~dT - “17“ ’ (8) 

where Zq = 0.015 is the Solar metallicity. Thus t z is the time 
scale on which the planet’s metal content increases by the 
amount contained in that planet at Solar metallicity. Since 


the fragment does not migrate it is not tidally disrupted in 
these isolated runs. 

Figure [6] shows the results of this series of runs for three 
initial planet masses, M p = 0.5, 1 and 2Mj. The top panel 
shows the mass of the core assembled inside the fragment by 
the time it collapses versus the metal loading time on the 
horizontal axis. For M p = 1 Mj, the core mass is essentially 
independent of t z . The bottom panel of Fig. [6] explains the 
reason for this finding to some degree. The panel shows the 
time it took the fragment to collapse. Unsurprisingly, the 
faster the metals are added to the planet, the faster it col¬ 
lapses, so that the core indeed has a shorter time window to 
grow. For M p = 1 Mj the two trends - the faster core growth 
and the shorter time available for that - nearly cancel each 
other out, yielding a nearly constant result. For M p = 2Mj, 
the dependence is mostly flat, but at the longest t z the core 
mass actually decreases. This is correlated with the time it 
takes for the fragment to collapse (the bottom panel of Fig. 
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EJ. Physically, for M v = 2 Mj, the longest t z planets collapse 
faster because their contraction is dominated by radiative 
cooling rather than pebble ac cretion. Radiative contraction 
is faster at lower metallicities dHelled fc Bodenheimerlfeoill : 
lNavakshinll2015ah . and hence the result. Finally, for the low¬ 
est mass planets in Fig. El the faster the metals are loaded 
into the planets (the shorter t z ), the lower the mass of the 
cores assembled. 

The different behaviour shown by the planets of differ¬ 
ent initial mass in their response to a variance in t z is due to 
a fact that a number of factors, rather than just one, control 
the core’s assembly rate inside the planet. The most impor¬ 
tant of these factors are: the central density of the planet, 
its metal abundance, gas temperature, and the core’s lumi¬ 
nosity as this affects the convective flux and the convective 
grain mixing which opposes grain sedimentation. 

The weak dependence of the core’s mass on the metal 
loading time (and hence on the metallicity of the disc), 
shown in Fig. El explains why the populations of the rocky 
planets in Fig. El are broadly distributed over metallicities. 
Indeed, if a roughly constant fraction of precollapse frag¬ 
ments behaving like those in the isolated simulations used 
for[6]is tidally disrupted, then there would be cores of similar 
masses at a range of metallicities. 

However, note that there is a precipitous fall in the num¬ 
ber of rocky planets at the highest metallicities in Fig.[o]that 
would not be explained in this picture. We believe that this 
decline in the number of rocky cores with 2 at high z has to 
do with a changing fraction of precollapse fragments expe¬ 
riencing tidal disruptions. The higher the metallicity of the 
disc, the less likely the fragment is to be disrupted. This can 
be seen from the bottom panel of Fig. El The fragment col¬ 
lapse times are quite short ( ;$ 10 5 years) at the low t z end. 
The fragment migration times are longer than that for most 
of our runs in the grid of models, and hence most of these 
fragments are expected to survive the tidal disruptions. 

This is why at the highest metallicities most of our pre¬ 
collapse fragments manage to collapse and either survive as 
a gas giant at the end of the simulation or reach the inner 
boundary of the disc. Either way, their rocky cores are never 
released from under the massive gas envelope, and this is 
why there are so few rocky planets at the highest metallici¬ 
ties in Fig. El This fall is probably too strong compared with 
observations. We come back to the issue in the Discussion 
section. 

6.2 Large separations 

Figure El shows the metallicity correlations for planets sur¬ 
vived the disc migration phase at the outer a > 10 AU region 
of the star. For the Jovian mass planets in the figure, the 
planets that are stranded far from the host star are those 
that migrate the slowest. Therefore, these fragments are not 
likely to be threatened by the tidal disruption, so can take 
as long as they need to contract. This is why the popula¬ 
tion of gas giants is now broadly distributed over the host 
star metallicity, with just a shallow minimum at around the 
Solar metallicity. This is quite different from Fig. El where a 
strong positive corelation with metallicity was present. 

The Super Earth and the Earth-like planets in Fig. El 
are the products of tidal disruptions of precollapse gas frag¬ 
ments, so these tend to come from those fragments that 



Figure 5. Number of gas giant, super-Earth and Earth-like 
planet groups versus host metallicity. Only planets left at sep¬ 
arations 0.1 < a < 5 AU are included in the figure. 


Tidal disruptions off 



10 4 10 5 

Metal loading time, t z [years] 


Figure 6. Top: the mass of the solid core assembled inside the 
precollapse gas fragment as a function of the metal loading time, 
t z , for non-migrating fragments of three different masses, as indi¬ 
cated in the legend. Bottom: The time it takes the fragments to 
contract from the initial state and collapse via H 2 dissociation, 
as a function of t z . 
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Figure 7. Number of gas giant, super-Earth and Earth-like 
planet groups versus host metallicity. Only planets found at sep¬ 
arations a > 10 AU are included in the figure. 

were migrating faster (physically ’’released to migrate” ear¬ 
lier, when their discs were more massive) than the gas giant 
planets in the figure. The metal correlation of the population 
of the Super-Earth planets in the outer 10 AU is not actu¬ 
ally that different from that in Fig. [5] with similar physics 
driving it. 

Earth-like planets are however much less abundant at 
high separations from the host, and especially at the highest 
metallicity bins. This is due to the fact that very few high 
metallicity fragments are disrupted in the outer disc. When 
they are disrupted, they tend to contain more massive cores 
because these cores had more time to grow by grain sedimen¬ 
tation than cores disrupted in the inner disc (which would 
usually migrate in much faster). This is another interesting 
example where predictions of TD do not follow the simple 
’’more metals more cores” logic. Hopefully this prediction 
may be observationally tested, although these relatively low 
mass planets are not likely to be found by direct imaging 
surveys in near future due to their exceedingly low lumi¬ 
nosities. 


7 DISCUSSION 

In this paper we argued for interpreting the metallicity corre¬ 
lations of observed exoplanets in the framework of the Tidal 
Downsizing hypothesis for planet formation rather than in 
the context of the Core Accretion model. In the latter the¬ 
ory, the fact that massive cores form at all metallicities but 
gas giants prefer metal-rich environments is odd. Once a 
massive core is assembled, accreting a massive gas envelope 
on the top of that core is actually easier in metal-poor envi¬ 
ronments. Lower opacities in the envelope allow the heat of 
the core’s assembly to escape from the gas envelope faster, 
accelerating its contraction and the eventual collapse. This 
can be seen from the well known result that the critical core 
mass for atmosphere collapse in the CA theory in creases ( al¬ 
though weakly) with opa city of the envelope le.g.. lStevensonl 
1 19821 : llkoma et al.ll2000ll . So, if massive cores are equally 
abundant at all metallicities, one would then expect more 


gas giant planets at low metallicities than at high metallic¬ 
ities, which is observationally not at all the case. 

In contrast, formation of gas giants in the TD hypoth¬ 
esis is not preconditioned by the assembly of solid cores. As 
we found here, there is actually no simple relationship be¬ 
tween these two types of planets in general. Both the rocky 
core dominated planets and the gas giants originate from 
the same source - the precollapse H 2 dominated fragments 
born by the gravitational instability of the outer massive 
disc - but the efficiency with which the different kinds of 
planets are assembled from this raw material is a function 
of a number of factors. 

Our numerical population synthesis models allowed us 
to study these issues in some detail in this paper. Focusing 
first on the planets surviving the disc migration phase in 
the inner 5 AU of the disc (SED. we found that increasing 
metallicity of the environment leads to a more ab undant for¬ 
mation of gas giant planets, as earlier found by iNavakshinl 
d2015bl l for coreless gas fragments. This result is driven en¬ 
tirely by the pebble accretion accelerating collapse of the 
fragments at high metallicites and helping more of them to 
survive. 

However, formation of solid cores in our models does 
not correlate in a monotonic way with the metal content of 
the parent discs for either Eartli-like or Super Earth plan¬ 
ets. There are two primary reasons for this outcome. Firstly, 
higher pebble accretion rates on a gas fragment do not nec¬ 
essarily make more massive cores. The core’s mass is equal 
to the integral of the core’s accretion rate over the duration 
of time period when the core can grow, e.g., / 0 ® row M c dt. 
The core’s accretion rate, M c , increases with increasing peb¬ 
ble accretion rate, but the duration of the core assembly 
phase, fgrow, decreases. For these reasons the core’s mass 
turns out to not depend sensitively on the pebble accretion 
rate (which is proxy for the disc metallicity), as could be 
naively expected (see Fig. [6]). 

Secondly, in accord with the explanation offered for the 
observed gas giant planet correlation with the metal abun¬ 
dance in TD, fewer precollapse fragments are disrupted by 
the stellar tidal forces at higher metallicities. Thus, not only 
the cores are not necessarily much more massive at higher 
metallicities, but they are also less likely to be released from 
inside the overlying massive gas envelopes. This is the rea¬ 
son why the population of the rocky cores nose-dives at the 
highest metallicity bin in Fig. [5] This particular result is 
however preliminary and may weaken if we are to consider 
earlier generation of gas fragments that may migrate more 
rapidly, so that a sizable fraction of precollapse gas frag¬ 
ments is tidally disrupted even at the highest metal abun¬ 
dances (which is not currently the case in our models). In 
that case there will probably be more rocky planets at the 
high metallicity end in Fig. [5] 

Finally, in H6.2l we compared the metallicity correlations 
of the simulated planets in the outer disc (a > 10 AU) to 
that found in the inner disc. Our models predict that gas 
giant planets of all metallicites may survive at these ’’cold” 
regions of the disc, since tidal disruptions of precollapse frag¬ 
ments are far less likely there. One interesting coincidence 
here is the fact that the best well known system of giant 
gas exopl anets probabl y formed by GI of the disc, HR 8799 
(e.g.. iMarois et al'1l2008l 'l . is metal poor ([Fe/H] « —0.5). 
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8 CONCLUSIONS 

We believe the observed metallicity correlations of plan¬ 
ets as a function of their size/radius give us a valuable 
clue as to how planets form. It is highly significant that 
gas giant planets are found almost always arou nd metal- 
rich stars (lGonzaledll999l : [Fischer fe Valentill2005l l but plan¬ 
ets smaller than ~ 4A® or less massive than Neptune are 
abundant around star s of any metallicity (ISousa et al .1120081: 
iBuchhave et alj|2012h . 

If solid cores are assembled first as a prerequisite to gas 
giant planet formation, then it is hard to see why there are 
massive cores in low metallicity environments but no gas 
giants since the lower the opacity of the envelope the easier 
it is to convert a massive core into a gas giant. 

On the other han d, if gas fragments are the nurseries 
of m assive solid cores llHelled fe Schubertl 120081 ; I Nayakshin] 
l201ll l. and if the fragmen ts are preferential ly destroyed in 
metal-poor environments (INavakshinl l2015aT ) . then the ob¬ 
served metallicity correlations are to be expected as long as 
the fragments are able to hatch massive cores before they 
are tidally disrupted. The population synthesis models pre¬ 
sented here appear to be consistent with the observations. 
They also predict that gas giant planets and Super Earths 
may be abundant at large separations from the star at all 
metallicities. 
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